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SUMMARY 


The plane wave propagation « the stability and the rectangular duct mode 
problems of a compressible inviscid linearly sheared parallel, but otherwise 
homogeneous flow, are shovm to be governed by Whittaker’s equation. The exact 
solutions for the perturbation quantities are essentially the Whittaker M-func- 
tlons where the non-dimensional quantities t, n and 4xn^ have 

precise physical meanings. A number of known results are obtained as limiting 
cases of our exact solutions. For the compressible finite thickness shear layer 
it is shown that no resonances and no critical angles exist for all Mach numbers, 
frequencies and shear lay.^r velocity profile slopes except in the singular case 
of the vortex sheet. 
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SYMBOLS 


a 

b 


c 


*>• 



f»g 


k.kjj 


m 


n 


-V 

n 


Speed of sound of the homogeneous fluid 
Velocity profile slope of the shear layer 
Disturbance phase speed (in general complex) 
Unit vector in the direction of the mean flow 
Independent solutions of equation (1) 

Wave numbers of incident wave 
Second index of the Whittaker M-funotions 
Index of refraction 
Wave normal unit vector 




w(r) 


x,y ,2 

z. 


Linear combinations of f and g in equation (5) 
Parameter of transformation in equation (1) 
z-component of the velocity perturbation 
Rectanular coordinates 
Shear layer thickness 


A,B,C,D 


A,B 


H 

K 


Wd 

t?f 

M(z) 


Functions of f, g and their derivatives in equation (7) 
Functions of A,B,C and D in equation (9) 

Heaviside function 

Inverse of the x-component of the disturbance Mach number 
Disturbance vector Mach number 

Mean flow vector Mach number 

Mean flow Mach number 

ll H ' 

Real part of 


Reflection coefficient 

^l'^2 phase component of R 

T Transmission coefficient 

Tj^,T 2 In and out of phase component of T 

U(z) Mean velocity 

W Dependent variable in Whittaker's equation 


V 


x-component of the wave vector of the pressure disturbance 
Non-dimensional varicible in equation (2) 

Angle of incidence of plane wave 

Independent variable in Whittaker's equation 
Non-dimensional parameter in equation (2) 

Perturbation velocity potentials 

Disturbance frequency 


INTRODUCTION AND BACKGROUND 

Studir.8 on comprosstble froa shear layers have not been extensive. In fact 
with tha exception of the earlier work of Graham and Graham (ref. 1) who studied 
sound propagation through a finite linearly sheared layer in the low-froguency 
limit, it has been only recently that Blumen et al (ref. 2) obtained an exact 
solution for tho stability of the shear layer with an hyperbolic tangent profile 
with the significant result that this shear layer is unstable to two dimensional 
disturbances for all Mach numbers whereas the vortex sheet is known to be un- 
stable only for M<2s/27 a result that cautions against modeling real sheared flows 
with vortex sheets, as has been the practice in a number of rocent noise re- 
saarch studies, since, as the authors point out, even the long wavelength char- 
cicteristics of finite thickness shear layers may be quite different from the 
corresponding properties of the analogous vortex sheet. In this study we con- 
sider sound propagation and stability in linearly sheared parallel compressible 
inviscid homogeneous flows. Work relating to the solutions of the pressure per- 
turbation equation has been that of Kuchemann (ref. 3) who also considered the 
stability of a boundary layer approximated by a linear velocity profile, the 
study of Pridmore-Brown (ref. 4) and that of Graham and Graham (ref. 1) . 

Kuchemann (ref. 3) obtained a formal series solution for the density per- 
turbation equation and he also arrived at a solution supposedly valid for large 
values of (our) parameter n *= i -M. His series solution, althoi 4 h it is given- 
in a cumbersome and lengthy form, is correct but his asymptotic solution is in 

serious error. Pridmore-Brown (ref. 4) solved the pressure perturbation equation 
in the short wavelength approximation, i.e., for large values of (our) parameter 

T = “ K. His asymptotic solutions may also be in error. Graham and Graham (ref. 

4b 

1) studied the problem of a plane wave incident on a linear velocity profile free 
shear compressible inviscid layer. They used entirely a series solution of the 
density perturbation equation which they independently rederived apparently un- 
aware of the earlier work of Kuchemann. Because they used the series solution 
only, they were unable to give proofs for the range of the parameters t and n 
for which ordinary, total or amplified reflection occurs, although correctly 
identified the regions intuitively. More importantly and for the same reasons, 
they could neither prove the existence or non-existence of reasonances or crit- 
ical angles, but for the case of "sufficiently thin — but not zero thickness — 
shear layer" nor could they draw any conclusions for either the large Mach 
number M or large x cases. 


*Work supported under NASA Grants NASA 2007 and NASA 676 to the Joint 
Institute of Aeronautics and Acoustics. 
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SOLUTION OP TUP PRESSURE PERTURBATION EQUATION 

In a homogeneous inviscid compressible parallel shear flow having a linear 
vvloctfcy profile in the z-direction only, i.e., U « U(z) » bz, it may be easily 
shown that, starting either from the linearized equations of motion or directly 
using the appropriate linearized form of the convective wave equation, the z- 
dependent part p(z) of the pressure perturbation p(^) is governed by the equa- 
tions 


^nn ■ i - i) P « o d) 

where 

n » ^ - M, 4 t * ^ K, and M « M(z) - ~ • (2) 

M is the local Mach number and K and u acquire the following meanings depending 
on the problem at hand: 

(i) Free Shear Layer ; Propagation of a plane wave of wave vector k and fre- 
quency 01 impinging on the shear layer from a half-space (z<0) of relative 

rest and at an angle 6 measured from the z-axis (- ~ < q < 4. i); 

2 — — 2 • 

K » sin9 

(ii) Free Shear Layer s Stability for assumed disturbances of the form 

p($)» p(z) eia(x-ct) ^ possibly complex); 

K » and 01 » Oc 


(iii) Sound Propagation in Rectangular Ducts t Modes for assumed disturbances 
of the fo™ p(ic) = p(z) eia(kx-wt) (jc and oi real) : 


1 

The transformation P “ w(5),C*qn with q = 4it, reduces equation (1) into 
Whittaker's equation for W, so that the two independent solutions f and g of 
equation (1) are; 


8 - 


(4iT)“'I“"' 





(3) 
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W*Mt© M at® the Whittaker M-functions and T,n»K are defined following equ- 
ation (2). It is easily shown that for t and n real^th® functions f and g are 
also real with f being even and g odd functions of n whereas both f and g are 
oven functions of t. Moreover series and asynptotic expansions for f and g are 
roadily obtainable from the known properties of the Whittaker M-functions (ref. 

5 and 6) . The series expansion agrees with the series obtained by Ku'chemann 
(ref. 3) and Graham and Graham (ref. 1) although the form that results from equ- 
ation (3) above is not only more compact, but also faster converging. The asym- 
ptotic forms of f and g and the^r derivatives with respect to n are obtained in 
terms of ^?m+l with m « ± using Oliver's method (ref, 5 and 6) and they 

are in disagreement with both the results of Kuchemann (ref. 3) and Pridmore- 
Brown (ref. 4). This was expected as mentioned in the introduction since 
Klichemann essentially seeking an expression for large n neglected 1 compared to 
rr* in the last term of our equation (1) which is tantamount to setting t « 0 in 

Wiittaker' s equation, whereas Pridmoro-Brown by applying banger’s method obtained 
only a non-uniform leading tern of an asymptotic expansion in terms of Airy fun- 
ctions. These and other details may be found in reference 7. 


THE FINITE THICKNESS: 

Plane VJave Propagation 

Wt* fonaider th« two-dimensional finite inviscid rompresslblo shear layer of 
thtekne^fi with velocity profile 

U •» 

“ bz 

>■ 0 0"*z 


in tlw (x,k) plane and a time-harmonic; monochromatic plane wave incident from 
the a 0 luilf-apace with wave vector 1* and wave number k * in the (rest) frame 
of reference of the stationary fluid at z<0, in an otherwise homotifeneous fluid 
in the entire (x,c) plane. The velocity potentials in the lower region {c.f re- 
lative rnst) and the upper region of uniform flow are,: 

. ^ I nk{xsinH+?.c:osO-at)_j^j^^-Uk{xslnfl-zcosO-at)l| ,z<p 


',|+j 


+ik 


xsin0±(z.-zj) |sin0| (n‘ 


■ 7 It 


z>z. 


O ** R.P.J+lTe" 
u 

where R.P. denotes "real part of"; the first term in equation (5) represents the 
incident wave coming f.vom the half-space z£0 and R and T are respectively the 
complex reflection and transmission coefficients for the velocity potential, and 
the upper signs are taken for n^O,* In the middle (shear layer) region with the 
velocity profile equation (4) the pressure perturbation p(r) and the z-component 
of the velocity perturbation are given by: 


w(r) 


p(r) * p^^^(n)sin |k(xsin0-at) 
a i..(l) 


w(r) 


'(n>- 


+ (n)cos 
1 

,< 2 ) 


k(xsin0-at) I - Pp (n)sin 


k(xsin0-at) 
k (xsin0 




where pf^)(n) and are linear combinations of the independent solution f and 

g, equation (3), of equation (1), i.o. 


(n) = 


^ 11 ^ 






+ 322 ^ 


(5) 


*This representation used by Miles (ref. 8) and Graham and Graham (ref. 1) is 
consistent with the radiation condition as postulated by Milos. Actually 
Sommerfeld’s radiation condition does not apply for plane waves and the diffi- 
culties arising in such a case have been discussed by Lighthill (ref. 9). At 
any rate these representations for the velocity potentials insure that the re- 
flected and transmitted waves are outgoing in a reference frame fixed in the 
upper fluid and are consistent with Milos’s postulate and Ribner's intuitive 

picture (ref. 10). 


4 


wHU j const^mtu. Writimi also R » ^ iR^# * • + iT., for thp cc'ri'U'X j.^*- 

{ li*. t i.»n anU trrnmiBelon roufficionfcs R anfl T rosppctively ati4 applyin^T the 
iHuinduiy s.*onvUti- ns {A’ontlimity of tho prc»f»ur« perturbatUm p(t'} anl - 

oni'nt wCr) of fcht* volov'tty pcrturbattonl at fehp Intorfaras z*0 and 

naparat iivq tonX and innginary partn In tho ronulting equations on> obtai*}‘i 
a uyntom of eight ilnuar algobrnic equations for the detormination of tV.o oigJ:r. 
unknowns! Rj^, After somewhat tediou* but straightforward algebra the 

t'ol lowing exprt*«s!ion« for the roflactlon and tranamiesion coefficients are ob- 
talmuls 


Tainaiff (aTb)^ + (cTd)^ 


tanO 


1 \/n 


I 


(« 


for lnj^|>l, - I <t»< Rj,, real, and 



1, 'f » Tj^ »■ Tg ** 0, 


for [nj^|<X» “ 


R^, Tj complex 


with, 


A >• 
B « 

C »» 


K 


^yi:n?r 


[■. 


(O)g^(l) - f^(l)g^(0) 


(4T)(sgnK>yn^Il f(i)g(0) - f(0)g{l) 


MIW I MI>_e»'W^ 


(sgnK)y(i 2 -i 


fU)g^(O) - f^(0)g(l) 


D » fj^U)g(o) - f(o)gj^(l) 


(7) 


with K ■ sinO 

The upper signs in equation (6) hold for n,>l and the lower signs for n <-l, 
in both cases jnpi. In equation (7) we used the notation 0 and 1 in the argu- 
ments vf f and g and their derivatives with the understanding that 0 designates 
evaluation at and 1 designates evaluation at n“r»i»nl _ « 

•L Z~ 2 «i| 


1 

sin*'* 


‘M 


i.G., at the two edges of the shear layer. 


It is clearly seen from equations (6) that the various reflection regimes 

ore ! 



t 

R^<1 

; Ordinary Reflection 

-i<tij< i-i 

f 

r2«i 

! Total Reflection 


$ 

R^>1 

: Amplified Reflection 
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I 

HieHO rogimast aro quit© .nnalogoua fco the thrae roglmai found by Kil#» (rof* 
B) 4Jid Ribnar (iref . 10) for tha vorfeax ahaefe caaa and Intufcivaly arrived at by 
Graham and Graham (rof, 1), The above conditions in aquation (B) imply that al“ 
though the values of the reflection and transmission coefficients depend on the 
frequency « and the velocity profile slope b, the dependence is only through the 
angle of incidence 0, and the two parameters t« ~ sinO and whereas 

«, tha conditions for the three reflection regimes are independent or w or b and de- 

pend only on which is the Mach number — of the x-component of the phase 

velocity of the incident wave front relative to the relative Mach number Mj^ of 
• the two uniform flows confining the shear layer. 


The limiting case of the vortex sheet is easily obtained in the limit t-K> 
(high-frequency or long acoustic wavelength limit) whereas the low-frequency or 
short wavelength limit is obtained by letting in equation (6); in the former 


case: 




K Hh4 



T*K) 


V'/hich agrees with the results of Miles (ref, f) and Ribner (ref. 10), whereas in 
the short wavelength case l-R^ becomes the Heaviside function H: 


1-r 2 * H (l-Hj ^) , T*>« 

Resonances and Critical Angles 

In this section we give a formal proof that in the amplified reflection re- 
gime, Hj^<-1, there are no resonances and in the ordinary reflection regime, 

there are no cricital angles. 

First it is easily deduced from equation (6) that excluding the singular 
cases of 0-K) or M-x» (in special ways) we may assume that neither A“C«0 or B*D*0, 
nor all four A,B,C and D may be zero simultaneously. It is next seen from equ- 
ation (6) that resonances exist if the denominator in the expression for is 
zero in the eunplifying regime i.e., if 

A = = 0, A» A+B and B *« C+D (9) 

with A,B,C, and D given by equations (7). But A in equation (9) above is just 
the deteeminant of the coefficients a^^ in equations (5) of the system of the 

four equations determining aj^j. Since only the first equation of that system is 
inhomogeneous with right hand side proportional to t, a solution for the a.^ 
exists if and only if either A « A^ 4- ^ O and t y 0 or A = 4 •+ O^and 

t-K). It thus follows that equation (9) has solutions, only for T’K), and this is 
precisely the limiting case of the vortex sheet; i.e., resonances are possible 
only for the vortex sheet. One may also obtain the same result by algebraic 
manipulation of the general expressions for R^^ and R^: 

2 2 2 2 

R a -C — , Rp «= Y (AD - BC) , A « (A 4 B)^ 4 (C 4 D)^ (10) 

■•■A A 


6 


for tht' case of amplified reflection (lower signs). 

For the critical angles we use the expressions in equation (10) with the 
upper signs (ordinary reflection n»>l). If critical angles exist, then R *R„“0 

X Xjjk 

and using equation (10) we may easily deduce that for Dj^O the ratio S’ « k 

may then only attain the value -1 for zero reflection. But this implies that 
A+B«C+D»0 which is precisely the condition for the existence of resonances equ- 
ation (9) which we have just shown that do not exist for a finite thickness 
shear layer. 


Tills riNlTE THICKNESS I.AVER 


Stability Consid«rationa 

For tho lay*r wquation (4) tha boundary valua ijroblam laada to tha following 
oguattons 1(A F B) i (C - D) • 0 ( 11 ) 

wlwro A,H,C, and U in gonaral oomplax ara givan by equation (7) with K 
v^ll’h («qn K) oimlttod in th» axpraaaions for H and C. The roots of the ibove 
oqiuttlan give the dopendenca of the phase speed c, or of the frequency \'n 
the wavenumber n. For temporal amplification a is real and positive, c atui 
u\«w,n'e real. For the neutral ntabllity line liowevur, in either case c»»Cj.tlc^ 

with c^**0, i.e., c is real, thus K is real and one may distinguish the cases 
K>1» In(>^ currosponding to supersonic (upper signs) or subsonic (lowor signs) 
disturbances! and relative Mach numbers respectively. 

Comparing equation (XI) and equation (0) we see that the resonances, were 
they to exist, would obey the system of equations A-fB»0, C+D*0 whereas the 
neutral stability characteristics are determined by the system of equations 
AFB»0, C-D«0, with A,B,C,D real. Thus, in general one does not 'xpect any con- 
nection between resonances and neutral stability eigenvalues except in n 
singular case of the compressible vortex sheet case which is discussed b&. 


Special Case: The Compressible Vortex Sheet 


As before, V!fo let t X) in equation (11) and (7) to obtain 

(1- 4»)f K2r,2 ^ (i-„?.)7 . 0 


( 12 ) 


vduno n and K may be complex. Excluding the singular rases K-»-l and nH, as well 
as K ■ -I* ►«', we consider the case where the square root terms in equation (7) 
have the name signs, i.o., the cases K<l,|nl>3. or K‘*l,ln|<l. Tho formal solution 
of equation ill 2) above is . 
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vv'hlch in in aqrnnwnnt wllh th<.» reaulfca of Umdau (raf. 11). It ia a amttor of 
jumplQ nlqobnn to i\\\t.i\4 fhnt in osdot to H.itlaCy the tuoqualltion jnl^l only 
tho \ippor (luinun) alun in tho aquaro rot>t turm in tho abovo ©qiuition^^Joi' thn 
i'iilonVitUios Hlumld bo lotAlmHl. 'rbnn two noutrnl oliionvalnos ato not por- 
wlnntblo which Is in-ooisely tha jroault of Milos 12) which ho avtivori at in 

a totally difforonfc way. namaly by conHldoring tho vottex sho©t ntabllltv as an 

initial value problem. It ia finally worth notinq tl»at for tho vortox shoot 
cane tho stabilrky equation (ll) for neutral eiqenvaluos, becomes C-b"0, wht'roas 
for the piano wave propagation case as we naw previously equation (*i) for the 
rosonancos becomes C+DK), since for the vortex sheet, and A and I) are 0(i) 

whereas C and D are 0(1). Thus it is only for the vortox sheet that resonances 

and neutral stability eignevaluas are given by the same equation i.o. *• D^. 
I’or the finite shear layer the roots of equation C*^)) and (11) the two equations 
aie in general difforont. Tn fact, we )»ave sliown that .ilthough equation (h) 
m.iy have real roots, t)\ere are no resonances Cor tlie finite thickness compress- 
ible shear layer. 


PHYSICAL MEANING OF , T and 4 th^ 

b 


Variable n » 

stability and rectui 
for n ; 


i -M. 


1/K is 


- ... for plane wave propagation, 

tgular duct mode studios respectively, and we may thus write 
■y .► 


"■d - ’V> 


■V 

.Of 


Thus n is a relative Mach nui'.be.* iieasure i v., it is the parallel to the 
$nean flow component of the disturbur.Ow- {phase speed based) vector Mach number 
relative to the (relative ) mean flow Mach number and thus it is a measure 
of the components of the relative speeds of the disturbance and the mean flow 
in the direction of the mean flow. 

Parameter ^ acquires the simple meaning of a characteristic Strouhal 
nvimber of the flow by writing; 

S »« (Shear layer thickness) x (Disturbance frequency) ^ zm „ tn 

2 Mean Plow Speed bz b 

Parameters t and 4xn : For propagation of a plane wave incident from a 

homogeneous half-space (z<o) at an angle 0^^ (-|l0ol+ |) z-axis, it is 

easy to show that the wave normals ii are independent of x, i.e., that all wave 
normals of a ^given^z-stratum are parallel. Thus defining an index of refrac- 
tion n “ ( 1 +M£*'u) a (1+Msin0)“*, it is easy to show that n « l-MsinG^ =» sin0oG» 
n stnsinOn. Using these relations we may write; 


Local disturbance wavelength _ 


Relative refraction index change 


Vnl/n 


(!) . n 

— sxn0„ 
b 0 


2 2 2 
n *=4x0, and 4xn is the 


argument of the Whittaker M- functions in our general solutions f and g of the 
pressure perturbation equation, t itself also attains the simple meaning 
1 m/b Characteristic Strouhal Number 


1/K 


Parallel component of tlie disturbance Mach number 


CONOIAU5IONS 


In thin i'apnv havo oxaminovl iuv.nv» nsponts of plant' wavn pi'opnsjat iv'n an.A 
'.lability of compronaiblo Invincid honKKii'nooun riowa c'haraot'C'i’l^.i'd by a I inoar 
v.'loaity profilo. Tho fooun nC thin Btiuly ban boon on tho noaroh for o.sai't 
nv'lul'tonr of the porturbalion oquationn which brinq foi'Hi tho naliont co;r,non 
toaturoa of all such parallel flows. The ossont ,il ronclur.ions of this stiKly 
are; (1) Tho a-clopemlont part of fcho prosauro (or domiity) dinturbanco of 
•nneh flows Is oovernod by Whittaker's equation witli independent solutions 

n'M (4irn^)» with m « 4 where M are tho Whittaker M-functions and n,:, 
i i , 'm 

and 4in‘ admit tho following interpretations: 

n « *’^f “ Relative Mach number parallel to mean flow 




tu/b 



Strouhal number /Disturbance Mach number component 
/ in the direction of the mean flow 


4ni'^ 


1 Vh\j/n 


Local 


d is turbanco v/.avo length 


/ Relative refractive 
index change 


C;l) solutions to a number of other parallel flow problems may be obtained as 
limiting cases from our exact solutions. Such flows include the compressible 
vortex sh(H't (t '0) , the incompressible vortex sheet n '0, ("n)>0), the in- 
compressible shear layer (i n) -finite) , and the short wavelength 

approximation of the compressible finite shear layer (r>«). (3) The compress- 
ihlo finite thickness layer has no resonances and no critical anglos for all 

Mach numbers, frequencies, shear layer thicknesses and shear profile slopes 
except for combinations of the singular values of 0 or w for, w and bj two 
such combinations (b*»«>, z, ►© but bz^ finite or constitute the compress- 

ible vortex sheet case. 
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